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RESUMEN 

Para la nebulosa planetaria NGC 7009, presentamos una comparacion de 
dos conjuntos de dates sin publicar del archive del Very Large Array temados 
con una separacion temporal de 8.09 anes para cenfirmar los movimientos prepios 
ebservados en el optico en sus ansas. Determinamos valeres de 23±6 y 34±10 mas 
afio~^ para las ansas este y oeste, respectivamente. Hay evidencia marginal de 
que las densidades de flujo de los chorros que conectan a las ansas con el cuerpe 
principal de la nebulosa disminuyeren en un 30% en el periodo que separa a las dos 
ebservacienes. Tambien establecemes un h'mite superior a la expansion del cuerpe 
principal de la nebulosa planetaria, obtenicnde un limite inferior de ~700 pc para 
la distancia a ella. 

ABSTRACT 

For the planetary nebula NGC 7009, we present a comparison of two unpub- 
lished Very Large Array archive data sets taken with a time separation of 8.09 years 
to confirm the proper motions of its ansae observed in the optical. We determine 
values of 23±6 and 34±10 mas yr^^ for the eastern and western ansae, respectively. 
There is marginal evidence suggesting that the flux densities of the jets that con- 
nect the ansae with the main body of the nebula diminished in about 30% over the 
period between the two observations. We also set an upper limit to the expansion 
of the main body of the planetary nebula, setting a lower limit of ~700 pc for its 
distance. 

Key Words: PLANETARY NEBULAE: INDIVIDUAL (NGC 7009) 
TECHNIQUES: INTERFEROMETRIC 



1. INTRODUCTION 

NGC 7009 (FN G037.7-34.5), the "Saturn Neb- 
ula" , is a well-studied elliptical PN, which has a jet- 
like system as well as a pair of low-ionization knots 
along its major axis. This type of knots, also known 
as ansae (AUer 1941), are now known to be present in 
many planetary nebulae and have received the names 
of bipolar, rotating, episodic jets BRETS (Lopez, 
Vazquez, & Rodriguez 1995; Lopez 1997), and fast, 
low ionization emission regions (FLIERS; Balick et 
al. 1993; 1994). There are several models proposed 
to explain the origin of these structures, but no one 
is generally accepted (for a discussion, see Steffen, 
Lopez, & Lim 2001). The radial velocity observa- 
tions of several of these objects by Balick, Preston, 
& Icke (1987) indicated velocities in the order of sev- 
eral tens of km s~^ with respect to the systemic ve- 



locity of the planetary nebula. On the other hand, 
proper motions of these knots have been measured 
only for a handful of objects (KjPn 8; Meaburn 1997, 
Hen 2-90; Sahai et al. 2002, NGC 7009; Fernandez, 
Monteiro, & Schwarz 2004). In particular, at radio 
wavelengths there are, to our knowlewdge, no mea- 
surements of proper motions of ansae in planetary 
nebulae. In this paper we present such measurement 
for the ansae in NGC 7009. We also discuss an ap- 
parent change in the flux density of the jets of this 
object and study the expansion of the main body of 
the nebula, setting a lower limit to its distance. 

2. OBSERVATIONS 

The two sets of observations used in this study 
were taken from the archive of the VLA of the 
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NRACQ- The epochs of the observations were 1989 
March 28 (epoch 1989.24) and 1997 April 29 (epoch 
1997.33), with a time separation of 8.09 years. Both 
sets of observations were made at 3.6 cm in the B 
configuration and each has a total of ^7 hours of 
on-source integration. In both epochs the source 
1331+305 was used as an absolute amphtude cal- 
ibrator (with an adopted flux density of 5.21 Jy) 
and the source 2131-121 was used as phase calibra- 
tor (with a bootstrapped flux density of 3.19±0.01 
Jy for the first epoch and of 2.68±0.02 Jy for the 
second epoch). 

The data were reduced using the standard VLA 
procedures in the software package Astronomical 
Imaging Processing System (AIPS) of NRAO and 
then cross-calibrated using the procedure of Masson 
(1986; 1989a; 1989b). 

3. INTERPRETATION AND RESULTS 
3.1. Proper Motions of the Ansae 

In Figure 1 we present 3.6 cm images for the two 
epochs analyzed. Three main structures are evident 
in the images. The first is the bright, main body 
of the nebula with an elliptical shape and angular 
dimensions of about 32" x 24". The second structure 
is constituted by the faint jets that emanate from the 
main body of the nebula and extend about 8" to the 
east and west. Finally, the third structure are the 
ansae, with angular dimensions of a few arc sec and 
that possibly are the termination points of the jets. 

The images were made with the ROBUST 
weighting parameter (Briggs 1995) of AIPS set to 
and a tapering of 150 kA in the (u,v) plane to smooth 
the angular resolution and enhance the signal-to- 
noise ratio of the relatively faint and extended ansae. 
The expanding proper motions of the ansae are evi- 
dent in the images. In Table 1 we present the values 
of these proper motions. The peak positions of the 
ansae have been obtained with the task MAXFIT of 
the AIPS software package. The errors in the po- 
sitions were estimated following Condon (1997) and 
Condon & Yin (2001). There are two previous re- 
ports of these proper motions in the literature. Filler 
(1965) estimated a value of '^16 mas yr~^ for the 
proper motion of the ansae in NGC 7009. A more 
modern determination is that obtained by Fernandez 
et al. (2004), who used HST images to obtain a value 
of 28±8 mas yr~^ for the eastern ansa (the western 
ansa was not included in one of the archive images 
they used). This value compares well with the value 

^The National Radio Astronomy Observatory is operated 
by Associated Universities Inc. under cooperative agreement 
with the National Science Foundation. 



of 23±6 mas yr~^ derived by us (see Table 1). As- 
suming a distance of 0.86±0.34 kpc, the average of 
14 available values (Acker et al. 1992; Fernandez et 
al. 2004), and taking into account the error from our 
proper motions as well as the relatively large error 
in the average distance, we obtain crude estimates 
for the velocities in the plane of the sky of 100±50 
and 140±70 km s~^, for the east and west ansae, 
respectively. 

If we assume that the ansae move ballistically (i. 
e. with no acceleration or deceleration) we estimate 
an age of ^^850 years if they originated from the cen- 
tral star and of ~300 years if they originated from 
the edge of the main body of the planetary nebula. 

3.2. A Decrease in the Emission from the Jets? 

There is marginal evidence in Figure 1 that, be- 
tween the two epochs observed, the emission from 
the jets that connect the ansae with the main body 
of the planetary nebula has decreased. For both the 
east and west jets, we find that the fiux density de- 
creased from ~0.9±0.1 mJy in 1989.24 to -0.6±0.1 
mJy in 1997.33, a decrease of about 30%. The ansae 
themselves, at the tip of the jets, do not show evi- 
dence of variation within the noise. 

The possible decreases in the fiux densities from 
the jets over the 8-year period arc marginally consis- 
tent with recombination theory. The electron den- 
sities of the jets have been estimated to be be- 
tween 1,000 and 4,000 cm~^ (Bohigas, Lopez, & 
Aguilar 1994; Bahck et al. 1994; Gongalves et al. 
2003; Sabbadin et al. 2004), implying recombina- 
tion timescales from 30 to 120 years. Evidence of 
variation in the surroundings of NGC 7009 has been 
found by Bohigas et al. (1994), who report a previ- 
ously undetected [SII] condensation about 20" to the 
NE of the main body of the nebula. On the other 
hand, as a whole, the planetary nebula has a fiux 
density of 654±3 mJy in 1989.24 and of 649±5 mJy 
in 1997.33, so the total flux density remains constant 
within <1 %. 

3.3. An Upper Limit to the Expansion of the Nebula 

Comparison of data taken with years or decades 
of separation often reveal the expansion signature 
in planetary nebula (e. g. Guzman, Gomez, & 
Rodriguez 2006). When subtracting the older data 
from the more recent one, the characteristic expan- 
sion signature shows as an outer region of positive 
emission, surrounding an inner region of "negative" 
emission. In Figure 2 we show an average image 
from the data of the two epochs (obtained by con- 
catenating both data sets) and a difference image, 



PROPER MOTIONS OF ANSAE IN NGC 7009 



3 



; 45 
50 
55 
22 00 
05 



: 








1 

1989.24 

: 






























1 



- -112130 



11.5 11.0 10.5 10.0 
RIGHT ASCENSION (J2000) 




11.5 11.0 10.5 10.0 09.5 
RIGHT ASCENSION (J2000) 



S 

a 43 



1 1 1 1 

1989.24 
Eastern Lobe 


II 














) M ^ 









II 



1997.33 
11 21 41 h Eastern Lobe 




21 04 12.60 12.55 12.50 12.45 12.40 12.35 12.30 12.25 12.20 
RIGHT ASCENSION (J2000) 



\ \ \ I \ \ I L 

21 04 12.60 12.55 12.50 12.45 12.40 12.35 12.30 12.25 12.20 
RIGHT ASCENSION (J2000) 



1989.24 
Western Lobe 




21 04 09.30 09.25 



9.20 09.15 09.10 09.05 09.00 08.95 08.90 
RIGHT ASCENSION (J2000) 



1997.33 
Western Lobe 




21 04 09.30 09.25 



9.20 09.15 09.10 09.05 09.00 08.95 08.90 
RIGHT ASCENSION (J2000) 



Fig. 1. 3.6 cm contour images of NGC 7009 for 1989.24 (loft) and 1997.33 (right). The top part of the panel shows the 
whole of the nebula, while the middle and bottom parts of the panel show the eastern and western lobes, respectively. 
The crosses mark the peak position of the lobes for 1989.24. Note the small displacement of the peak position of the 
lobes for 1997.33. The contours are -4, 4, 5, 6, 8, 10, 12, 15, 20, 30, 40, 60, 80, 120, 200, 300, 400, and 500 times 19.5 
/ttJy beam^^, the average rms noise of the two images. The synthesized beams (l!'74 x l!'50 with a position angle of 
-8° for 1989.24 and l!'76 x l!'38 with a position angle of -15° for 1997.33) are shown in the bottom left corner of the 
images at the top of the panel. 



4 



RODRIGUEZ & GOMEZ 




-11 21 30 

35 



o 
o 
o 

CM 
-3 



g 



40 
45 
50 



uj 55 



22 00 
05 

21 04 12.5 



1997.33-1989.24 







12.0 11.5 11.0 10.5 10.0 
RIGHT ASCENSION (J2000) 



09.5 09.0 



Fig. 2. 3.6 cm contour images of NGC 7009 for all the data (top) and the difference 1997.33-1989.24 (bottom). The 
contours are -4, -3, 3, 4, 5, 6, 8, 10, 12, 15, 20, 30, 40, 60, 80, 120, 200, and 300 times 12 /xJy beam"^ for the top image 
and -5, -4, -3, 3, 4, and 5 times 20 /iJy bcam^^ for the bottom image. The synthesized beam (l"00 x 0"81 with a 
position angle of —10°) is shown in the bottom left corner of the images. 
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TABLE 1 

PROPER MOTIONS" OF THE ANSAE 



Ansa 






'^totalO 


P.A. (°) 


Id (mas yr ^) 


Eastern 


+0.14±0.G5 


+0.11±0.05 


0.18±0.05 


52±16 


23±6 


Western 


-0.23±0.08 


-0.14±0.08 


0.27±0.08 


239±17 


34±10 



'^AnA is the displacement in right ascension, Aoec is the displacement in declination, Atotal is the total 
displacement, P.A. is the position angle of the displacement, and fi is the total proper motion. 



1997.33-1989.24 (obtained by subtracting the clean 
components of the 1997.33 image from the (u,v) data 
of 1989.24). Those images were made with the RO- 
BUST weighting parameter of AIPS set to and 
no tapering. Although there is evidence of resid- 
ual emission in the brightest parts of the nebula, we 
failed to detect a clear expansion signature (see bot- 
tom of Fig. 2). An upper limit to the expansion can 
be obtained following Gomez, Rodriguez, & Moran 
(1993) and Guzman et al. (2006). We assume that 
the nebula expands in a self-similar way by a factor 
of (1 + e), where e << 1. From the non detection 
of an expansion pattern, we estimate that e < 0.007. 
A larger value of e would produce residuals in the 
difference image larger than those observed. A lower 
limit to the distance, D, of the planetary nebula is 
then given by 



D 

pc 



> 0.211 



km 



At 



yr 



where is the average angular radius of maximum 
emission, Vexp is the expansion velocity of the nebula 
at the radius of maximum emission, and At is the 
time interval between observations. NGC 7009 is 
an elliptical nebula and we approximate the angular 
radius of max;imum emission as the geometric mean 
of the major and minor semiaxes. From Figure 2, 
we estimate 6 ~ 12". Prom Weedman (1968) we find 
that at this angular distance v,,.j.p ~ 35 km, s^^. We 
then obtain D > 700 pc, a lower limit consistent with 
the weighted average of 14 available values (Acker et 
al. 1992; Fernandez et al. 2004), 860±340 pc. It 
should be noted that the expansion parallax distance 
technique can be affected by the motion of the shock 
or the ionization fronts (Mellema 2004) and that our 
lower limit could be underestimated by 20-30%. 

4. CONCLUSIONS 

We analyzed archive VLA observations of the 
planetary nebula NGC 7009 for two epochs sepa- 
rated by 8.09 years. Our main conclusions are: 



1. We measured the proper motions of the ansae, 
obtaining values of 23±6 and 34±10 mas yr~^ for the 
eastern and western ansae, respectively. This is the 
first time that proper motions of ansae in planetary 
nebulae are determined with radio observations. 

2. There is marginal evidence suggesting that the 
fiux densities of the jets that connect the ansae with 
the main body of the nebula diminished in about 
30% over the 8.09 years spanned by the observations. 

3. We failed to detect the expansion of the main 
body of the planetary nebula, obtaining a lower limit 
of ~700 pc for its distance. 
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